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. . . A measurement of the relative branching fraction of Bl? — !> J/?/) /o (980), /o (980) —>■ tv^-k~ to 

I I B° ^ J/ip(t>A -> K+K- is presented. The decay mode B^l J/i/'/o (980) is an interesting mode 

. since it is a CP-odd eigenstate which could be used in CP-violating studies. Using approximately 

' 8 fb~^ of data recorded with the DO detector at the Fermilab Tevatron Collider, a relative branching 

CN \ fraction of 0.210 ± 0.032 (stat) ± 0.036 (syst) is found. 

Oh' 

C/^ \ I. INTRODUCTION 

m : 

^-H , The CP-violating phase in mixing mixing, has been measured [l|, [4I using J/4'4' decays. The 

measured absolute value is larger than predicted by the Standard Model (SM) Q, but is statistically consistent 

^ with it. The decay products in 5° J/if'fo{980) are in a CP-odd eigenstate and can provide a more direct 

' measurement of this CP-violating phase. Measuring this CP-violating phase using — > J/'0/o(98O) decays 

JL , mode can aid in reducing its uncertainty. 

Based on estimates the relative branching fraction should be large. Using hadronic Df decays, Stone and 

(~| ' Zhang [1, Q estimated the relative width to be: 

- ■- ^ „ r(i3g ^ J/^/o(980);/o(980) ^ tt+tt-) ^ 

>■ ^= r(Ro^ j/j,^.^^ K+K-\ 



The LHCb collaboration has reported [oj a first measurement of i? = 0.252j;o o32to o33. The Belle col- 
laboration has made a measurement of the branching fraction B{B'^ — ^ J/'0/o(98O); /o(980) — tt+tt~) = 
• (1.16to;?9(stat.)to;i7(syst.)to;i8(Ng(.)g(.))) x IQ-^ [13]. The CDF collaboration has also measured the rel- 
. ative branching fraction and finds /?=0.257 ± 0.020 (stat) ± 0.014 (syst) ll]. This article provides a new 
' measurement of the relative branching fraction using the DO detector collecting data at the Fermilab Tevatron 
I . Colhder. 

^ ' This note provides a new measurement of the relative branching fraction from DO. 
• ^ . 

X ■ 

^ ■ II. RELATIVE BRANCHING FRACTION 

! 

To determine an absolute branching fraction, various efficiencies, branching fractions, and cross sections need 
to be known, as well as the integrated luminosity. However, by measuring a relative branching fraction, several 
terms common to both the B'^ — > J/-0/o(98O) branching fraction and the i?° — )■ J/iIj(I) branching fraction cancel 
giving: 

„ ^ g(i?g ^ J/^/o(980); /o(980) ^ ^+Tr-) ^ iVBg^j/^/o(980) x efl^^^^^ 

BiBO ^ J/^0; ^ K+K-) Nso^j,^^ X ■ ^ ' 

All that is required to measure a relative branching fraction are the relative yields and the relative reconstruction 
efficiencies of the two decay modes, e^i^'^^^'^ and e^i^'' ^^^^^^"^^^^ . 



III. SELECTION CUTS 



A. Analysis Cuts 

The data set of an integrated luminosity of approximately 8 fb~^ was divided into four periods corresponding 
to different detector configurations called Runlla, Runllbl, RunIIb2 and RunIIb3. 

The initial sample of i?^ — J/V'/o(980) was found by first reconstructing J jif) — >■ [i^ ^J.~ candidates by requiring 
that two oppositely charged muon candidates with transverse momentum > 1.5 GeV form a common vertex. 
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Since the DO detector has a hmited abihty to separate kaons from pions, ah reconstructed tracks not associated 
to a. J/ip are considered for reconstructing /o(980) and (p candidates. The tracks are assigned the pion mass 
when searching for B'^ J/ipfo{980) and the kaon mass when searching for B'^ — ?> J/ipip. Two tracks with a 
minimum pr of 300 MeV, having an invariant mass 0.7 GeV < M^+^- < 1.2 GeV, and being consistent with 
coming from a common vertex were considered as /o(980) candidates. Finahy, the /^+/i~7r+7r~ candidates were 
required to have a common vertex and have an invariant mass between 5.0-5.8 GeV. 

Similar requirements were apphed to the initial sample of B'^ — >■ J/ipc/) candidates. The only different re- 
quirements were that 0.91 GeV < Mj^+j^^ < 1.05 GeV and the fj,^iJ,~K^K~ candidates were required to have 
an invariant mass between 5.0-5.8 GeV. Due to the invariant mass requirements on M^+^~ and Alj^+j^-, two 
tracks cannot be considered both a /o(980) and a candidate. The final data sample was then formed by 
applying the additional requirements: 

• All runs without optimal performance of muon, silicon microstrip and central fiber trackers are omitted . 

• All events that only fired a trigger that required muons with a large impact parameter were removed. 
J ftp selection: 

• Both muons are required to be detected as a track segment in either one or three layers of the muon 
system and be matched to a central track. 

• At least one muon must be detected as a track segment in three layers of the muon system. 

• Both muons must have at least one hit in the silicon microstrip tracker. 

• 2.9 GeV < M^+^- < 3.2 GeV 
/o(980) (0) selection: 

• Both pions (kaons) from the /o(980) {(p) candidate must have at least 2 hits in the central fiber tracker. 

• Both pions (kaons) from the /o(980) (0) candidate must have at least 2 hits in the silicon microstrip 
tracker. 

• Both pions (kaons) from the /o (980) (0) candidate must have at least 8 hits total in the silicon microstrip 
tracker and the central fiber tracker. 

• The momentum of the leading pion (kaon) from the /o(980) (0) candidate must be greater than 1.4 GeV. 

• /o(980) {(p) candidate pt must be greater than 1.6 GeV. 
selection: 

• 0.91 GeV < M^+„- < 1.05 GeV (when searching for J/V'/o(980).) 

• 1.01 GeV < Mk+k- < 1-03 GeV (when searching for J/^Jcj).) 

• pt(SO) > 5.0 GeV 

• Proper decay length [l3|, L, significance, L/cf{L) > 5, where (j{L) is the uncertainty on the proper decay 
length. 

B. Boosted Decision Trees 

It is known that boosted decision trees (BDT) [13, [ll] are a powerful tool for separating signal from back- 
ground. Signal and background samples are used to train the BDT and a discriminant is determined for each 
event. By making a selection on the value of the BDT discriminant, the signal to background ratio can be 
vastly improved. We use the Monte Carlo (MC) PYTHIA program jl6| to generate B'^ and the evtgen program 
[l7j to simulate its decay. Two MC background samples were produced: a prompt sample (directly produced 
J/ip) and an inclusive sample (all decay processes B*^ — J/^ + X). A MC signal sample of i3° — >■ J/'0/o(98O) 
events was then used to train the BDT on both the prompt and inclusive background. A BDT discriminant 
was found for both the prompt and inclusive sample and used in the analysis. A total of 36 different kinematic 
variables were used to train the BDT consisting of isolation variables, transverse momentum of the daughters 
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FIG. 1: BDT distribution after training for both signal (blue) and inclusive background (red). 



and grand-daughters of the and vertex quality of the i?" and its daughters. Figures [T] and [5] show the BDT 
distributions for the training and test samples for the inclusive and prompt background. 

The BDT cuts were determined only using the 1 fb~^ of Runlla data. A narrow window around the nominal 
/o(980) mass was chosen to keep the signal to noise ratio high. Using a mass cut of 0.96-1.0 GeV on the tt+tt" 
mass, the BDT cut value was chosen where both S/a/B and the signal yield were high. In this way, the BDT 
discriminant for both the inclusive and prompt BDT was required to be greater than 0.35. 



IV. YIELD RESULTS 



A clear peak is found when the tt+tt" invariant mass is near the nominal /o(980) mass. It is expected 
that the B^ signal can be fitted to a Gaussian distribution, which provides a fitted mean mass (/j,) and width 
(<t) for the B^ peak. Since backgrounds are large, a cut of ±2tT around the fitted B^ peak is used to identify 
the /o(980) mass peak. A clear /o(980) mass peak is observed when the /i+/i~7r+7r~ invariant mass is within 
±2(7 of the fitted B^ mass, see Fig. [31 To decide on a tt+tt" mass window to use for this analysis, a fit to the 
/o(980) mass peak is performed. The /o(980) has a large width [3l and is just under the KK mass threshold. 
This changes the line shape from a simple Breit Wigner form, particularly for higher masses and so the tt+tt" 
mass distribution is fitted using a functional form based on Flatte [19J, convoluted with a Gaussian function, 
that takes into account the opening of the KK threshold. The lineshape found from fitting the /o(980) in MC 
is used to fit the data. A tt+tt^ invariant mass cut of 0.91-1.05 GeV is applied to identify B° — 7> J/-0/o(98O) 
and is shown in Fig. |4l The B^^ J /ipf o{980) mass distribution was fit to a Gaussian signal with a background 
function consisting of a second-degree polynomial and a Gaussian function at lower invariant mass to take into 
account partially reconstructed B decays. 

Using identical cuts (except for the cut on the mass), a clear J/ip(j) peak is found and is shown in Fig. [5] 
Since the (f> peak is so narrow, the backgrounds are much smaller for — )■ J/ip(p. 

An unbinned likelihood fit was used to determine the candidate yields in each sample. The fit to the 
J/?/'/o(980) mass distribution shown in Fig. |4] gives the following results (statistical uncertainties only): 

B° mass = 5.3747 ± 0.0036 GeV; a = 0.0290 ± 0.0044 GeV; 498 ± 74 ^ 7/^-/0(980) candidates. 

The K^ K~ mass distribution was fit for a i?^ — > J/ip4' signal using a double Gaussian function with a 

second-order polynomial background. A fit to the J/4'(j) distribution shown in Fig. [5] gives the following results 
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FIG. 2: BDT distribution after training for both signal (blue) and prompt background (red). 
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FIG. 3: 71"*"^" invariant mass distribution peaking at the /o(980) mass when the J/ip-K^-K~ mass is ±2a" around the 
fitted B1 mass. 

(statistical uncertainties only): 

mass = 5.3631 ± 0.0008 GeV; 2863 ± 61 ^ J/^cf) candidates. 

V. EFFICIENCIES 

To determine the efficiencies of the analysis, MC signal samples were used. To take into account the effects 
of the instantaneous luminosity, the MC samples were overlaid with zero bias data collected during each run 
period. In the generation of both the J/4'4' and the J/'0/o(98O) signal MC's, a preselection requirement of 
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FIG. 4: /i /i TT TT mass distribution peaking at the Bl mass when the yr vr mass is between 0.91 and 1.05 GeV 
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TABLE I: The reconstruction efficiency for B° -5- Jl%l}(j) and B° J/V'/o(980) for various running periods. 



Sample 


total reconstruction efficiency 


B° J/i()(j} Runlla 
B° J/TP(j) Runllbl 
B° J/TP(f> RunIIb2 
B° J/V'/o(980) Runlla 
B° J/tljfoi980) Runllbl 
B° -> J /ipf 0(980) RunIIb2 


0.0231 ± 0.0004 
0.0191 ± 0.0004 

0.00636 ± 0.00018 
0.0191 ± 0.0004 
0.0146 ± 0.0003 

0.00529 ± 0.00015 



TABLE IL Reconstruction efficiencies for different run periods. 



Run period 


Relative reconstruction efficiency , , 

^ J/ ijfQ{9SQ) 
^TFrn 


Runlfa 
Runllbl 
Runllb2 


1.21 ± 0.03 
1.31 ± 0.04 
1.20 ± 0.05 



Pt > 0.4 GeV was demanded on both kaons (pions) from the 0(/o(98O)). Since the pr distributions for the 
pions and kaons may be different, the preselection efficiencies of this cut must be determined. To determine the 
preselection cut efficiencies, two additional MC sets were also generated with no pT cuts on the pions (kaons). 
By comparing these two results, the preselection cut efficiencies were determined. 

We found that the reconstruction efficiencies depended heavily on the MC sample used since the instantaneous 
luminosity was different for the various run periods, therefore we determined the reconstruction efficiencies for 
each run range separately. The instantaneous luminosities for runs taken during RunllbS were similar to the 
instantaneous luminosities for runs taken during RunIIb2 so the reconstruction efficiencies found from RunIIb2 
were used for RunllbS. Table U shows the results on the efficiency analysis using MC signal samples. Table 
|T] shows that the absolute reconstruction efficiencies vary in each run period, however Table show that the 
relative reconstruction efficiencies are relatively stable. However, the differences in the relative reconstruction 
efficiency is considered a systematic uncertainty on R. 

VI. SYSTEMATIC UNCERTAINTY STUDIES 

background studies 

One possible peaking background that affects the — > J/'0/o(98O) yield measurement is the non-resonant 
— > J / ijjTT'^ 'iT~ background. This background was studied by measuring the yields in 7r"'"7r~ invariant mass 
less than the /o(980) mass. The tt+tt" mass distribution from 5° J/V'tt+tt" background where the tt+tt" 
are non-resonant should have a much broader distribution, so determining the yield for lower tt+tt^ masses 
will allow a determination of the contamination in the /o(980) signal region. 

In determining the tt+tt^ mass window to study, it is important to choose a window where one does not 
expect other resonances (i.e., — > J/ipK*). The tt+tt^ mass window of 0.8-0.9 GeV was chosen since in this 
mass range there should not be any B^ — > J/t/jK* events. In fitting the distribution for any possible signal, the 
signal /i and a are constrained to be the values found from the fit to the 5° mass in the /o(980) signal region. 
The fit yields 80 ± 75 events, giving no statistically significant evidence of any B'^ J/i/^tt+tt^ non-resonant 
background, so no correction was applied. 

B. Analysis cut variation 

To cross check that the results do not vary with the exact value of the analysis cuts, the choice for each 
analysis cut was varied around its nominal value. This is an important test since the selection criteria was 
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TABLE III: Fractional change due to varying the exact choice of analysis cuts on the relative branching fraction 



Cut 


e {J/tp4>) 




event yield —>■ J/ip4> 


event yield B'^ —J- J/ipfo 


effect on R 


BDT inc > 0.3 


1.000 


1.017 


1.020 


0.958 


0.96 


BDT inc > 0.4 


0.993 


0.980 


0.975 


0.945 


0.98 


BDT pro > 0.3 


1.000 


1.002 


1.000 


1.007 


1.01 


BDT pro > 0.4 


1.002 


1.000 


1.000 


0.991 


0.99 


Pt{B1) > 4.5 GeV 


0.997 


1.000 


1.000 


1.000 


0.99 


Pt(B°) > 5.5 GeV 


1.000 


0.995 


1.000 


0.952 


0.95 


Pt(/o(980)) > 1.0 GeV 


1.000 


1.000 


1.000 


1.000 


1.00 


Pt(/o(980)) > 2.0 GeV 


1.000 


0.987 


1.000 


0.980 


0.99 


tt/K pt > 1.0 GeV 


1.210 


1.099 


1.172 


1.133 


1.06 


Tv/K PT > 1.8 GeV 


0.724 


0.771 


0.797 


0.744 


0.88 


L/a{L) > 4 


1.057 


1.047 


1.056 


1.035 


1.01 


L/a{L) > 6 


0.946 


0.951 


0.944 


0.967 


1.02 



TABLE IV: Effects of changing the fitting choices 



Parameter 




J/i/'/o(980) yield 


Nominal fit (Gaussian signal + second order polynomial background with fit range 5.1-5.8 GeV) 




498 ± 74 


Third degree polynomial background 




446 ± 72 


Background function exponential+polynomial 




423 ± 67 


Fit range 5.1-5.6 GeV 




437 ± 78 


Fit range 5.15-5.8 GeV 




427 ± 63 


Fit range 5.05-5.8 GeV 




449 ± 71 



determined with 1 fb^^ data from Runlla, and it is important to verify that this did not introduce a bias into 
the measurement. Table Hill shows the results from this study. As can been seen from the table, the value of R 
does not depend significantly on the exact choice of selection requirement. 

C. Fitting cross checks 

Due to large backgrounds arising from combinatorics and partially reconstructed B decays, there are signifi- 
cant uncertainties in the exact background shape. Therefore different parameterizations were used to describe 
the background and different fit regions were used to fit the data. The background polynomial was changed 
from a second-degree polynomial to a third-degree polynomial. The fit range was changed from the nominal 
5.1-5.8 GeV and finally a different functional form for the background was used by changing the background 
shape to a polynomial plus an exponential. 

As can be seen from Table IIVI there is a fairly large variation in the number of signal events for _B° — >■ 
J/'!/'/o(980), indicating that the background shape is difficult to model. This fitting systematic gives the largest 
systematic uncertainty on R. A study was performed using same-sign pions and forming the mass distribution 
from fi^ fj,~ TT^ . However, it was found the the same sign pion distribution did not describe the measured 
background and so could not be used to help constrain the background shape. A similar study of varying the 
fitting choices was performed on the — >■ J/ip'P sample, however since the backgrounds are much smaller and 
easier to describe the event yield numbers changed by less than 1%. 

A summary of the uncertainties on the BR are summarized in Table |Vl 
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TABLE V: Statistical and systematic uncertainties in branching fraction ratio, R 



Source 


Uncertainty 


Statistical 
Systematic from fitting 
Systematic from different MC samples 


0.f49 
O.fSO 
0.0858 



VII. FINAL BRANCHING FRACTION RATIO 



The decay 5° — >■ J/tpfo{980) is an interesting decay mode since it can allow a measurement of the CP-violating 
phase in mixing. 

A measurement of the relative branching fraction using approximately 8 fb^^ of data yields: 



R 



B{B"^ ^ J/V./o(980); /o(980) ^ t^+tt- 
B{BO ^ J/tl}(t)- 4) K+K-) 



0.210 ± 0.032 (stat) ± 0.036 (syst). 



The relative branching fraction of 5° J/i/;/o(980), /o(980) tt+tt" to 5° J/iPcf), (jj K+K' should be 
large enough to allow a measurement of the CP-viiolating phase in 5° mixin using the decay 5° — J/'i/'/o(980). 
An analysis to measure (j)s using the decay B° — >■ J/'(/)/o(980) is currently being pursued. 
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